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We develop further our extension of the Ellis-Bruni covariant and gauge-invariant formalism 
' to the general relativistic treatment of density perturbations in the presence of cosmological 

magnetic fields. We present detailed analysis of the kinematical and dynamical behaviour of 
■ perturbed magnetized FRW cosmologies containing fluid with non-zero pressure. We study the 

C*~) , magnetohydrodynamical effects on the growth of density irregularities during the radiation era. 

Solutions are found for the evolution of density inhomogeneities on small and large scales in the 
presence of pressure, and some new physical effects are identified. 
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i . 1 Introduction 

5h . 
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In a recent article, |IJ], we examined the behaviour of cosmological density perturbations in a uni- 
verse containing a large-scale primordial magnetic field, by means of the Ellis and Bruni covariant 
and gauge- invariant approach Q. Our assumptions were: first, that the conductivity of the medium 
is infinite; and second, that the background universe, though permeated by a coherent magnetic 
field, remains spatially isotropic to leading order. The first approximation is a standard simplifi- 
cation of Maxwell's equations, ignoring any large-scale electric field while preserving the desired 
coupling between matter and the magnetic field. The second approximation was introduced at a 
later stage of our analysis, to allow the direct comparison between our results and those from previ- 
ous Newtonian treatments. Starting from a general, inhomogeneous and anisotropic, cosmological 
model we provided the exact, fully non-linear, evolution formulae for all the basic gauge- invariant 
variables. These equations are valid irrespective of the field's strength and can be linearized about 
a "variety" of smooth background universes. In Q, the main objective was to establish a fully rela- 
tivistic treatment of magnetized density perturbations. To achieve this, we focussed upon the dust 
era and compared our results to those obtained earlier by the Newtonian treatments of Ruzmaikina- 
Ruzmaikin, tj, and Wasserman, 0. In addition to the desired agreement with the non-relativistic 
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analysis, our method suggested weak corrections to the evolution of density disturbances. These 
are generated by both the isotropic and the anisotropic pressure that the magnetic field introduces 
into the cosmological model. Our results confirmed the relative unimportance of the field for the 
evolution of super horizon-sized density disturbances. 

The present article extents our study to the radiation era by including the isotropic pressure of a 
perfect fluid in the calculations. During this period the kinematic evolution of the universe follows 
more complicated patterns than in the dust era. The sources of the extra complexities are changes 
in the fluid motion due to its non-vanishing pressure, which manifest themselves in a number of 
ways. For instance, the acceleration of the fluid depends now on pressure gradients as well as on 
the spatial variations of the magnetic field, with the two of them not necessarily acting in the same 
sense. So, unlike the pressure-free case, the field vector and the fluid acceleration are generally not 
orthogonal. Moreover, the time derivative of the magnetic field is no longer a spacelike 3-vector. 
All these factors leave their trace on the kinematics of the universe and further complicate its 
dynamic evolution. It should be emphasized that we do not consider plasma-physics complexities 
induced by dissipative stresses in the energy momentum tensor of the fluid. A separate study of 
this magnetohydrodynamical problem, both at linear and non-linear level, is possible due to the 
conformal invariance of the stress tensor and can be found in Subramanian and Barrow B . These 
studies provide a mathematical and physical basis for the evaluation of the observational effects of 
cosmological magnetohydrodynamics. The existence of magnetic fluctuations can leave observable 
traces on the structure of the microwave background on small angular scales, which might be 
detected by future satellite missions if the magnetic field is strong enough to influence the formation 
of large scale structure. Subramanian and Barrow || have shown that, for a tangled magnetic field 
with present strength Hq ~ 3 x 10~ 9 G, one can expect an RMS microwave background anisotropy 
signal of order or larger, depending on the angular scale. The anisotropy in hot or cold spots 
could be several times larger. The formalism we have developed can also be used to trace the effects 
of the damping of magnetic field fluctuations on the photon and neutrino spectra emerging from 
the radiation era. On superhorizon scales the field evolves as though quasi-homogeneous and will 
create small expansion anisotropics which will produce anisotropies in the microwave temperature 
distribution. The existing observational data allows us to place upper limits of6.8xlO- 9 (ft ^ 2 ) 1/2 G 
on the present strength of any uniform (i.e. spatially homogeneous) component of the magnetic 
field §, if]. 

In this report we examine the case of a perturbed spatially flat Friedmann-Robertson-Walker 
(FFRW) magnetized universe filled with a single barotropic perfect fluid and derive the linear 
equations that determine its evolution. We consider the evolution of the basic kinematic and 
dynamic quantities and the magnetohydrodynamical effects upon them. We show that, in regions 
of subhorizon size, gradients in the energy density of the field, together with those in the fluid, 
decelerate the universal expansion and act as sources of positive spatial curvature. It also appears 
that the field tends to smooth out the curvature of the underlying 3-surfaces and, depending on 
the spatial curvature, it can act as a source of slightly accelerated expansion. 

Given the current interest in structure formation, we provide a set of four linear first order 
differential equations that governs the linear evolution of A, the scalar variable that determines 
the gravitational clumping of matter and describes the formation of structure directly. During the 

1 Q.o is the present value of the density parameter and h is the Hubble constant in units of lOOfcra/ sMpc. 
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radiation epoch the energy density of the fluid and that of the magnetic field fall as the inverse 
fourth power of the scale factor. As a result, the Alfven velocity is time-independent and our 
system of four differential equations becomes easier to handle. We provide analytic solutions at 
both limits of the wavelength spectrum and compare them to those of a non-magnetic universe. 
We show that although the density inhomogeneities retain their basic evolutionary patterns, the 
role of the field as an agent opposing their growth is clear. In particular, large-scale perturbations 
undergo a power-law evolution, similar to that of the non-magnetized case, but their growth rate 
is reduced by an amount proportional to the field's relative strength. At the opposite end of the 
spectrum the density contrast continues to oscillate. Here, the extra magnetic pressure has simply 
reduced the oscillation period. The conclusion is that during the radiation era the magnetic effects 
are just supplementary to those induced by the pressure of the relativistic fluid. After the radiation 
era ends, the field becomes the sole source of pressure in regions exceeding the associated at the 
time Jeans length. We find, in agreement with our earlier conclusions (see Q), that any large- 
scale magnetic influence ceases completely by the later stages of the dust era, although at earlier 
times it could have forced the density contrast to enter a brief oscillatory phase. The negative role 
of the field is also confirmed in the subhorizon regions. On scales between the Jeans length and 
the horizon, the magnetic field slows the power-law growth of the inhomogeneities by an amount 
depending on its relative strength. 

In [Q] we considered a general FRW universe, allowing for spatially open and closed unperturbed 
backgrounds. However, it is important to recognize that the gauge invariance of the magnetic field 
perturbations holds if and only if the underlying spatial sections are flat. Accordingly, all equations 
in must be linearized about a FFRW universe and every variable representing spatial curvature 
should be treated as a perturbation. As a result, terms in the linearized formulae of (l| that contain 
the background 3-curvature constant are non-linear and can be dropped at first order. This does 
not effect the results presented there but restricts their validity to almost-FFRW models. Notice 
that the gauge invariance of the magnetic field gradients still holds within a perturbed Bianchi-I 
universe due to the latter's spatial flatness. In Appendix A we give a full account of this question. 

2 Preliminaries 

2.1 Kinematic Variables 

Following H and []!(]] , we assume that the average motion of matter in the universe defines a future 
directed velocity 4-vector, Uj, corresponding to a fundamental observer (uiU 1 = — c 2 ), and generates 
a unique splitting of spacetime into "time" and "space" (1+3 decomposition). For any tensorial 
quantity T, the directional derivative T = T-^u 1 = ii*VjT denotes differentiation along the fluid-flow 
lines. The second order symmetric tensor hij = gij + UiUj/c 2 projects orthogonal to U{ onto what is 
known as the observer's instantaneous 3-D rest space We also introduce ®Vi, the covariant 

derivative operator orthogonal to U{ (^Vj/ijfc = 0), by totally projecting the corresponding 4-D 
operator. This is not, however, a derivative on a hypersurface unless the fluid flow is irrotational. 
Nevertheless, we will call the 3-gradient "spatial" for simplicity. 

2 At every event along the worldline of a fundamental observer, Sj_ is the normal to Ui 3-D subspace of the 4-D 
space tangent to that event. 
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The kinematic variables are established by decomposing the covariant derivative of Ui into its 
spatial and temporal parts. In particular, we have 

VjUi = + Uij + ^hij - -^aiUj, (1) 

where = ^Vqu^ — @h{j/3 is the shear {pijU 1 = OijV? = 0, = 0), ujij = ^^VuUji is the 
vorticity (uJijU 1 = uJijU J = 0), O = Vjii* is the volume expansion, and a% = Uj = u 3 V jUi is the 
acceleration (a^it* = 0). The magnitudes of the shear and the vorticity are a 2 = a^a^ /2 and 
to 2 = LUijLU^ /2 respectively. The expansion scalar, 0, defines a representative length scale (5) along 
the fluid flow by means of S/S = 0/3. In a non-rotating universe (i.e. when Uij = 0), U{ is a 
hypersurface orthogonal field and £j_ becomes a 3-surface, namely the instantaneous rest space of 
all the fundamental observers. 



2.2 Spacetime Geometry 

The global geometry of the spacetime is determined by the Riemann curvature tensor, conveniently 
expressed by the decomposition 

Rijkq — Cijkq 2 {9ikRjq ~i~ 9jq^ik 9jkRiq QiqRjk) g (.9ik9jq 9iq9jk) Ri (2) 

where Cijkq is the Weyl conformal curvature tensor, Rij = R^ ]k is the Ricci tensor and R = R l i 
is the Ricci scalar. Both the Ricci tensor and the Ricci scalar are determined locally by matter 
through the Einstein field equations. Conversely, the Weyl tensor describes long range gravitational 
effects, such as those of tidal forces and gravitational waves. By definition Cijkq satisfies all the 
symmetries of the Riemann tensor and is also trace-free. It decomposes into an electric and a 
magnetic part according to, pi , 



Cijkq ~ c2 (9vsr9kq P t V^srVkqpt) ME* ^ ( V ^ S r9kq P t + QijsrVkqpt) vfvPH*, (3) 

where 

9ijkq = 9ik9jq 9iq9jki (4) 

rjijkq is the totally antisymmetric spacetime permutation tensor and Ey, are respectively the 
electric and the magnetic components of the Weyl tensor. The latter have nothing to do with actual 
electric or magnetic fields but derive their name from the Maxwell-like equations they comply with 



[12]. Also notice that 



C ijkq = # |^ = ° o (5) 
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2.3 The Electromagnetic Field 

The electromagnetic field is represented by the antisymmetric Maxwell tensor F{j. This splits into 



an electric and a magnetic 4- vector respectively defined by, [10], 



Ei = F ijU J, and H { = ^-r, ijkq u^ F k ^ . (6) 

The above confirm that Eiu 1 = Hiu 1 = 0, which in turn mean that both fields lie on E 2 = E^E 1 
and H 2 = HiH 1 respectively denote the magnitudes of the electric and the magnetic fields. 

2.4 The Material Component 

As in Q, we consider a universe filled with a single perfect fluid of infinite conductivity (i.e. 
a = cEiJ l /E 2 — > oo, with Jj representing the current density). We can now drop the electric field 
from Maxwell's equations, which reduce to0 

2u i H i = ec, (7) 
^ Uj (a k H q -c 2 V q H k ) = c 2 K r P, (8) 

^Hi-^Hi = 0, (9) 
cr 

29) ■ \ ■ ■ • 

u\+u\-^-h l AW = h\W, (10) 



3 3 3 V 3 

where u>i = f]ij kq u J uj kq /2c is the vorticity vector and e = — JiU l /c 2 is the charge density. For our 
purposes the last two equations are the important ones. More specifically, (||) provides the familiar 
vanishing 3-divergence law for the magnetic field (i.e. ^ViH 1 = 0), whereas (jToj) , when contracted 
with the magnetic field vector, gives aijH % H 3 = 2QH 2 /3 + (H 2 )'/2 and simplifies the expression 
for the energy density conservation law, O. 

The energy-momentum tensor for a magnetized single perfect fluid of infinite conductivity has 
the form, [EL 




Tij = ( A 4 + 7j~2 ) UiUj + { p + — ] hij + Uij (11) 



with the pressure (p) and the mass density (/z), the latter including contributions from the internal 
thermal energy, related by a suitable equation of state. The symmetric, traceless and completely 
spacelike tensor Hy = H 2 hij/3 — H{Hj describes the anisotropic pressure induced by the magnetic 
field.0 

3 the merit of the infinite conductivity assumption is that, based on Ohm's law, it can accommodate a zero electric 
field with non- vanishing spatial currents (i.e. h? Jj 7^ 0). The latter condition is essential if one wishes to preserve 
the coupling between matter and magnetic field (see Appendix B in JlJ). 

4 In |l[ we represented the anisotropic magnetic stresses by Mj 3 - instead of 11^ . Other changes relative to that 
article are; at has replaced in as the acceleration vector; the 3-Ricci scalar has changed from K. into K, while now 
K, = S 2 K; the gradient ^ViH 2 is represented by Hi and not by Bi as in |lj; and the scalar B is no longer the 
Laplacian (3) V 2 if 2 but equals the dimensionless ratio 5* 2(3) V 2 H 2 /H 2 . 
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2.5 Inhomogeneity Variables 

In a FFRW universe all physical quantities are functions of cosmic-time only, while the shear, the 
vorticity, the acceleration, all the anisotropic stresses, the curvature of the spatial sections and the 
electric and magnetic components of the Weyl tensor vanish. Within a nearly-FFRW universe, 
spatial perturbations in the energy density and the pressure of the fluid, in the expansion and in 
the magnetic field are described by four key variables covariantly defined in Q and pjj. These 
are: the comoving fractional orthogonal spatial gradient of the energy density, Di = S^'V ifi/ n; 
the orthogonal spatial gradient of the pressure, Yi = n^Vip, where k = 8irG/c 4: is the Einstein 
gravitational constant; the comoving orthogonal spatial gradient of the expansion, Z = S^VjS; 
and the comoving orthogonal spatial gradient of the magnetic field, M.^ = kS^V jHi, with M.? = 
0, as (H) requires. Each one of these 3-gradients vanishes in a perfect FFRW model (see also 
Appendix A), thus satisfying the criterion for gauge-invar iance, |l3|]. Three additional gauge- 
invariant variables, which play a crucial role in our analysis, are the divergence of the acceleration 
(i.e. A = Via 1 ), its spatial gradient (i.e. Ai = ^ViA), and the spatial gradient of the curvature 
scalar K associated with Sj_ (i.e. Ki = ^Vj-KT). 

It is convenient to introduce the following local decomposition for the spatial gradient of Di, 

A 4j = S^VjDi = W i:j + S y + ±Ahi, , (12) 

where Wij = Ar^-i contains information about the rotational behaviour of Di, E»j = — A^hij/3 
describes the formation of anisotropies (e.g. pancakes or cigar- like structures), and A = Aj l is 
related to the spherically symmetric gravitational clumping of matter. Although in a general 
perturbation pattern we expect turbulence (i.e. Wij ^ 0) and anisotropic structures (i.e. Sjj / 
as well as material aggregation (i.e. A > 0), it is the latter scalar which is crucial for the structure 
formation purposes. 

3 The Linear Regime 

In reality [/,, p, Q and Hi must have a spatial dependence as well as a temporal one. Moreover, 
K, ai, dj, uJij, Hij, Eij and Hij will generally take non-zero values. Assuming that the observed 
universe is close to a FFRW spacetime, we can linearize the evolution equations by treating all the 
gauge-invariant quantities, along with their derivatives, as first-order variables. The exact, fully 
non- linear formulae have already been derived in [Q]. Here we give their linearized versions only. 

3.1 Evolution Equations 

The linear regime is monitored through the following combination of propagation formulae and 
constraint equations: 

(i) The conservation laws of the energy and the momentum densities of the fluid, respectively 
expressed by 

^ + (1 + ^)6 = 0, (13) 
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and 

K.ll.(~\ -i- una; 4-V: — 

s 



ku(1 + w)ai + Y % - \M m W = 0, (14) 



where the ratio w = p/uc 2 evolves according to 

w = -(l + w) I ^-w\ G, (15) 




with c 2 = p/fi representing the adiabatic sound speed. Although generally w is allowed to vary, 
when it remains constant along the fluid-flow lines (i.e. when w = 0) equation (|l!]) suggests that 
w = c 2 s /c 2 = constant, provided of course that 0/0. 

(ii) The propagation equations that determine the kinematics of the universe. These are Ray- 
chaudhuri's formula, 

B + % + ^- (1 + 3w) - A - Ac 2 = 0, (16) 

where A is the cosmological constant, and the propagation formulae of the vorticityQ and the shear 
tensors, respectively given by 

20 

Uij + -g-^ij = ^VyOi], (17) 

and 

29 f os A, kc 2 ^ 2 ^ / n 

where = Cikj q u k u q /c 2 is the "electric" part of the Weyl tensor. To first order, the scalar A is 
given by the 3-divergence of the fluid acceleration (i.e A = ®V*Oi). 

(iii) The linear propagation equations of D{, and M.%j, respectively governing the growth of 
spatial inhomogeneities in the energy density of the fluid, 



A = wQDi -(l + w)Zi- ——jM , ; // ' + 4 en, (19) 



2G . . „ A 2SQH 2 

in the expansion scalar |] 



~-Zi - ^y-A - 3c 2 A%]IP - c 2 ^^^' + 5^, (20) 



5 Equation (L7|) monitors the model's rotational behaviour through the vorticity tensor, as opposed to the vorticity 
vector used in [y (see eqn (89) there). Recalling that Ui = r/ijkqU-'uj q /2c the equivalence of the two formulae becomes 
evident. 

6 In El , based on the weakness of the magnetic field, we ignored the linear effects on the evolution of the expansion 
and the 3-curvature gradients resulting from the field's contribution to the active gravitational mass of the universe. 
Here we fully incorporate these effects via the second last terms in the right hand side of Jit) ) and (E?]) (compare 
them to eqns (91) and (99) in Notice that these quantities provide all the coupling between the magnetic and 
the matter inhomogeneities that is left, once the infinite conductivity approximation is abandoned in favour of a 
pure source- free magnetic field (see appendix B in Q). Although they make no qualitative difference and introduce 
negligible quantitative changes, both terms are included here for completeness. 
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and in the magnetic field vector, 

Mij = -^M ij -^H i Z j + KSH m y j (a ik + u:i k )- l ^ r {2H l a j +aiH j ) + 



3c 2 



a k H K hij + KhfR kqjs H q u s , (21) 



recalling that A{ = ^ViA by definition. 

(iv) The evolution of the magnetic field is governed by the four decomposed Maxwell's equations 
(see eqns ([7|)-([To|)), of which only 

ViH 1 = \a t H\ (22) 



hPH j = ( a% + oA - ^-h% ) W, (23) 



and 

are crucial for our analysis. The former verifies that the magnetic field is a "solenoidal" (i.e. 
^ViH 1 = 0), and the latter, when contracted with the field vector, provides a radiation-like linear 
evolution law for the magnetic energy density 

H 2 = ^, (24) 

where EI = 0. 

(v) We close this section with a brief discussion on the geometry of Sj_, the observer's instan- 
taneous rest space. Its curvature is characterized by the scalar 

K = 2 ^/ic 2 - ^ + , (25) 

so that 2 /3 = ft/uc 4 + Ac 2 to zero-order. When there is no vorticity and only then, K coincides 
with the 3-Ricci scalar of the spacelike hypersurfaces that define the instantaneous rest space of all 
the fundamental observers. Its propagation formula, 

suggests that in the linear regime the fluid acceleration acts as the sole source of spatial curvature 
through its 3-divergence. Following [O], we describe the spatial variations of the 3-curvature by 
the gauge-invariant vector C% = S 3 Ki and provide a supplementary relation between Di, Zi and 

48 9 2 

Ci = 2 Kf ic 2 S 2 Di + 2S 2 M ji W - ~^~Z h (27) 



which, by means of (|i9|)-(|2l|), leads to 



46 S 3 

Q = — g-A 4 . (28) 
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The above propagation formula is consistent with equation 
the interdependence between the spatial curvature and the acceleration of the fluid flow. The 
advantages of choosing Cj, instead of K+, to describe the spatial variations in the 3-curvature, will 
become clear later. 



26) and, together with (H2), confirms 



4 The Case of a Barotropic Perfect Fluid 

Among the propagation formulae given above, which refer to a general perfect fluid with pressure, 
there is no equation for the evolution of 3-gradients in the pressure. The reason is that the prop- 
agation of Yi will be determined directly from (|l9|), once the material content of the universe has 
been specified. 

4.1 Equation of State 

Here, we extend the analysis presented in by considering a universe filled with a single barotropic 
perfect fluid. Its equation of state is[| 

p = p(jj), (29) 

suggesting that Vr^Vji// = 0. Consequently, the relation between pressure and energy density 
gradients becomes 

SYi = Kfic 2 s Di, (30) 
since (? s = dp /dp, relative to the observer's rest frame. 

4.2 Kinematic Evolution 
4.2.1 The Acceleration 

The energy density conservation law of the barotropic fluid is still expressed by 
momentum density conservation law, (H), together with (pC|), gives 

ai = n T" \ c 1 (— M M Hj ~ c l D i) > 

(1 + W)S \K/J, 1 J / 

for the acceleration of a fundamental observer. It depends both on gradients in the energy density 
of the fluid and on gradients in the magnetic field. Thus, the geodesic flow can still be preserved 
provided that the field gradients counterbalance those of the material component. The necessary 
and sufficient condition for this to occur isf] Di = 2M. uj-i W / k\m? s . 

Unlike the pressure-free case (see [|J), the acceleration of the barotropic fluid is not always 
normal to the magnetic-field vector. Alternatively, one might say that, when p / 0, the time 

7 From now on, all our results will refer to a barotropic fluid unless otherwise stated. We will also ignore the 
entropy contribution to the fluid- pressure. 

8 The geodesic flow condition can simplify the evolutionary relations of section 3.2 considerably. However, it does 
not appear to be consistent and we will not pursue the matter any further here. 



13). However, the 
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derivative Hi no longer lies on Ej_. We verify these statements by simply contracting ( |31| ) with Hi. 
We find that 

OjiT = = - f* WDi, (32) 
(1 + w)S 

where generally HiD % ^ 0. Obviously, on and ilj remain orthogonal if HiD 1 = 0. We can modify 
this condition by setting E = H 2 / /ic 2 , taking its 3- gradient and then contracting with the magnetic 
field vector. The result, 

WDi = —W^ViE, (33) 

suggests that the acceleration of the fluid flow remains normal to the magnetic field if and only if 
the directional derivative H l ^ViE vanishes (i.e. when the energy density ratio E does not change 
along the magnetic field lines). In this case, the time derivative of the magnetic field lies on the 
observer's instantaneous rest space. This fact can simplify, among other, calculations involving 
commutations between the spatial gradients of Hi. 

Spatial gradients in the fluid acceleration affect the expansion dynamics directly (see eqns 



Ql6|), (|17|) and (18)), as well as the spatial geometry (see eqn (Pq )). Consequently the following 
new decomposition of the acceleration's 3-gradient is of major importance. It is obtained directly 
from equation ([H]) via the commutation laws for the 3-gradients of scalars and spacelike vectors 



(see eqns (101) and (102) in Appendix B), relations (|i2[), (P3[) and the relativistic expression 



Uij = H 2 hij/3 - HiHj for the magnetic anisotropic stresses (see §7.4.1 in @ for more details). 

c 2 ( A \ AQH 2 H 2 

(3)V ^ = - (irt)ffl & + w - + 3 h V ~ 9^(1 + w f ij + 3Rtt^) (3) ^ - 

1 -( 3 )V/ 3 )V^ 2 + 1 - H k <$V k Mi3, (34) 



2[i(l + w) J Kfi(l + w)S 

where ( 3 ) Rij is the 3-Ricci tensor of the spacelike regions (given by eqn (83) in 0). In what follows, 
the trace, the skew part and the symmetric part of the above will be employed to analyze the 
magnetohydrodynamical effects upon the kinematics and the spatial geometry of our cosmological 
model. 

4.2.2 The Deceleration Parameter 

To begin with, the trace of (^ 

A = ( 3 )VX = - ^ A + - f K - - H " — B, (35) 
{l + w)S^ 3/x(l +w) 2fi(l + w)b^ 

where B = S 2 ^V 2 H 2 /H 2 , is substituted into (|16|) to produce Raychaudhuri's equation for a 
magnetized universe filled with a single barotropic perfect fluid of infinite conductivity. This formula 
is recast into the following alternative expression 

- q = -£-(! + 3w) - - - K + JT —^ U A + ^B) - A, (36) 
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where q = —SS/S 2 is the "deceleration parameter" (note that A = S 2 ^ V 2 ^/^ to first order). 
Clearly, the sign of the quantity on the right hand side of (|36| ) determines whether the expansion 
slows down or continues unimpeded. Not surprisingly, a spherically symmetric increase in the 
energy density of the field (i.e. B > 0), together with any material aggregation (i.e. A > 0), slows 
the expansion down. Their combined effect is of the first order and, as the Laplacians verify, it is 
confined to regions well within the horizon. However, while the energy density of ordinary matter 
(i.e. w > —1/3) always adds a positive value to the deceleration parameter, the contribution of the 



magnetic energy density depends on the geometry of According to (36), the coupling between 
the field and the 3-curvature slows down the expansion of spatially open almost-FRW universe (i.e. 
when K < 0) but accelerates perturbed Friedmannian cosmologies with positive spatial curvature 
(i.e. K > 0). This rather unconventional magnetic effect is global, though still first order in 
magnitude since K = in the background. It vanishes when the perturbed universe retains its 
spatial flatness. 

4.2.3 The Vorticity Tensor 

According to (|l7|), only the antisymmetric part of ( pi] ) affects the vorticity propagation. Since 
Wij = -(1 + w)QS 2 uj ij /c 2 to the first order (see @) and ^ Vp® Vj]H 2 = AQH 2 ^ /3c 2 , as the 
commutator of the 3-gradients of scalars (see eqn ( |101| ) in Appendix B) and the last of Maxwell's 
equations (see eqn (|23|) ) imply, we obtain 



of) / "if2 i | 

*« + T 1 - 2# h = ^(T7^ ,3,v ^i- (37) 



Notice that a cosmic magnetic field influences the vorticity of the universe solely through the 
antisymmetric part of the gradient field Mij, which itself describes the rotational behaviour of the 
magnetic field vector (see Appendix C.l in (lj). Also, according to p?]), the field has no effect at all 
when the directional derivative H k ^Vk-M[ij] vanishes, that is when curLf/j does not change along 
the magnetic field lines. 

4.2.4 The Shear Tensor 



The symmetric part of (|34|) together with its trace allows us to recast equation (|l8[), for the linear 
evolution of the shear tensor, into the following form 

^ + - -ir^^-MiW( ,3,v '- 3>v «-^ <3,v2 ^ 2 + 

h 2 f(^ n \ kc 2 

R (ij) ~ ~o hi J I + — Hi,- + 



^^^V k M {ij) -^ Ei , (38) 

Clearly, the shear anisotropies evolve in a rather complicated way under the simultaneous influence 
of a number of sources. According to ( |38| ) such sources are: the fluid; the magnetic field; the 
geometry of the observer's 3-D rest space; and the long range source-free gravitational field. The 
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magnetic influence is multi faceted. In particular, anisotropic spatial variations in the energy density 
of the field have a similar effect to those in the energy density of the fluid (the latter represented 
by Sjj). Also, the magnetic energy density couples with anisotropies in the spatial curvature to 
create an additional effect. Notice that any anisotropic patterns in the distribution of the magnetic 
field vector (described by M-Uj) since Ai? = 0) exert no influence at all if the directional derivative 
H k ^V 'fc-M(y) vanishes. 

4.2.5 The 3-curvature Scalar 

In the linear regime, the curvature scalar of the observer's instantaneous 3-D rest space evolves 
according to equation (pq) . Substituting the trace of (|3^ ) into the latter we obtain 

(39) 





or, since H 2 //ic 2 (l + w) <C 1, 

(40) 

Therefore, on regions of subhorizon size, any spherically symmetric spatial increase in the energy 
density of the fluid (i.e A > 0), or of the magnetic field (i.e. B > 0), acts as a source of positive 
curvature. This is a first order effect similar to the small-scale magnetohydrodynamical impact 
on the expansion of the universe (see the last term in the right hand side of eqn (|36|)). The 
magnetic influence also results in a global effect of the opposite type. As ( |39| ) reveals, the field 
tends to smooth out the curvature of the spacelike regions through its coupling with the background 
expansion. This "magnetic smoothing" , which here is of negligible magnitude, is analogous to the 
field's global magneto-geometrical impact upon the universal expansion (compare to the second 
term in the right hand side of (^6|)) both qualitatively and quantitatively. 

We close this section with a sort comment on the non-local magnetic effects illustrated in 
equations (|36|) and (|^). We attribute such behaviour to the vectorial nature of the field, as opposed 
to the scalar nature of quantities such as the energy density of the fluid or its isotropic pressure. 
Being a vector, the field interacts with the curvature of the spacelike regions (e.g. through the 3- 
Ricci identity) and this dependence creates the aforementioned effects. However, we do not suggest 
that any perturbed spacelike vector would have a similar impact. The magnetic influence outlined 
above depends crucially upon the specific properties of the field, as these are reflected in Maxwell's 
equations, and also in the unique way general relativity describes the magnetic anisotropic stresses 
(see comments on the derivation of eqn (|34])). 



4.3 Dynamic Evolution 

4.3.1 The Growth of the Inhomogeneities 

The introduction of a barotropic fluid modifies the key linear propagation equations ( |I9| ) , ( |20| ) , 
and (|2§|), More specifically, by means of (|3l"l), the former becomes 



2G 

Di = wQDi - (1 + w)Zi 2 M M R3 - ( 41 ) 

K/jJiC 
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The direct barotropic influence on the evolution of the expansion gradients comes through the 
spatial gradient of A,f\ 

l ~ (l + w)S V l + 3^1 + w)S 3 1 2fi(l + w) 1 c 2 Vj ^' [ ] 

where Tii = ^ViH 2 . Substituting the above into and using (|27|) we obtain 



Z t = -^ Zl -^D l -3c 2 M M W-c 2 M Jl W--^^V 2 D l 



s .(3) V 2 W ._2^©5;(3) v .^- ; (43) 



2/x(l + iu) C 2 

By means of (Q) and the fact that h? RjkU k = (see eqn (63) in the last term in (^), 
which describes the effects of spacetime curvature upon the evolution of magnetic inhomogeneities, 
becomes 

KhfRkqjsHOu 3 = KhfC kqja Wu s . (44) 

Moreover, using decomposition (||) of the Weyl tensor we may recast the above into 

KhfCw&u' = -Kr/ ik qs H k u q H S] , (45) 

showing that only spacetime ripples caused by the long-range gravitational forces, here represented 
by the magnetic part = i q i ^ <1 C kq j s vPu s /2c 2 of the Weyl tensor, affect the propagation of spatial 
inhomogeneities in the cosmic magnetic field. The above result together with equations (31) and 
|) allows us to transform (|2l"|) into 



J^+ w) ( 2H i D i + H J Di ~ H k D k h tj ) + KS Vi « s H k Uq H sj . (46) 

According to (|4l] ) it is the only the contraction AdujiHi that contributes to the linear growth of 
spatial inhomogeneities in the energy density of the medium. Therefore, taking the skew part of 
([46[) and then contracting with the magnetic field vector, we obtain^ 

2B 2k 

M m Hi = -—M {ij] W - —H^W + KSh { k h« (a ks + uj ks ). q H S W + 

K f s@ HuDftW, (47) 
3c 2 (l + w) [l 11 ' y J 

which will prove useful later. Notice the absence of the spacetime curvature term from the right 
hand side of equation (^7j). As a result, long range gravitational forces have no linear effect on the 
evolution of magnetized density perturbations. 

9 In deriving (^ij) we have treated ' 3 - ) ViW and ' 3 - ) ViC^ as first order gauge-invariant quantities. Though the former 
is straightforward to prove, the latter requires the gauge-independence of ' 3 'Vip and ^Vifi to be shown first. 
10 See §5.4 and §7.5 in jlf| for a detailed derivation of the complete set of the exact and linear propagation equations. 
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As far as spatial inhomogeneities in the curvature scalar are concerned, equations (|30|), (j3l|) 
and ( |4^ ) reshape their propagation formula, (28), into the following 



.(3)v 2 Dj + 



26S 3 



.(3) V 2 ^ + 



8£e 2 s 3 



3c 2 (l + u;) v ~* ' 3/xc 2 (l + ^) v "* ' 3c 4 ' 
implying that the gradient field Cj is invariant on large scales (i.e. when the Laplacian terms 



(48) 



are negligible) if ®V 
Maxwell's equations become, 



0. Finally, under the barotropic fluid assumption, reflected in (p2j), 



c 2 (l +w)S 



H L Dj 



and 



Hi 



(Jij + 10 ij 



2G 



hij ) W 



9/1 i \ o 3 



(49) 



(50) 



with the non-zero right hand side of (49) and the last term of 
in the fluid motion relative to the pressureless case. 



being direct results of the changes 



4.3.2 The Growth of the Density Gradient 

The dynamics of the inhomogeneity variable Di is governed by ([II]) , together with equations ([l3|)P 
and (|47]), or by the linear second order differential equation, which follows from (|4l|), by means of 
@, ll), dH), dH), © and (H). This equation is 

„2 
„2 

5w Ac 2 Di + 




— Mm4\IP + 



flC z 



(51) 



This is the generalization of formula (26) in [14] for a magnetized almost-FFRW universe. It 
has the form of a wave equation with extra terms due to the universal expansion, gravity, the 
cosmological constant, the magnetic field and the vorticity. The difference in the vorticity terms 
between equation (|5l|) above, and its corresponding formula (115) in ||, is due to the residual 
coupling between the divergence of the vorticity tensor and the energy density of the field that 
remains when p = 0. 

11 Alternatively, one can use instead of (^), on substituting Zi by d in (^lj) from (^]). 
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5 The Scalar Variables 



So far we have considered the evolution of gauge-invariant vector variables and in particular the 
propagation of Di, the gradient field that describes orthogonal to the fluid flow variations of the 
energy density. However, regarding the growth (or decay) of density inhomogeneities, the vector 
field Di contains more information than actually required. We can extract the information we need 
by adopting the local decomposition (|l2|) . Of the three additional variables mentioned there, the 
scalar A = S^V l Di (alternatively A = (S* 2 //i)( 3 ^V 2 /i to first order) is the most important one 
when addressing the problem of structure formation. 



5.1 Definitions 

Focusing upon A, which describes spherically symmetric spatial variations in the energy density of 
the matter, we also consider the following complementary scalar variables, 

Z = S^V*Zi, B = 4i (3) V 2 # 2 , (52) 

respectively related to spatial gradients in the expansion and the energy density of the magnetic 
field, and 

K = S 2 K, (53) 

representing perturbations in the spatial curvature. Notice that all but Z are dimensionless vari- 
ables. Also, B describes spherically symmetric spatial variations in the energy density of the 
magnetic field and it will be treated as the magnetic analogue of A. 



5.2 Evolutionary Equations 



The propagation equations associated with the above defined scalars (see §7.6 in [15] for details on 
their derivation) are 



A = wQA - (1 + w)Z 



<~>H 2 

3/ic 2 



B. 



(54) 



29 „ k/ic 4 . 
-Z — A 



3 

KC 2 H 2 



1 + w 



-(3)V 2 A 



KC 2 H 2 



B 



H 2 



2/x(l + w 



_(3) V 2 B 



(55) 



B 



4c 2 9 
3c 2 (1 + w) 



A 



-Z 



AQH 2 
9fic 2 (l + w 



-JC, 



and 



K 



4c 2 G 



-A + 



2QH 2 



3c 2 (1 + w) 3^c 2 (l + w 



-B. 



(56) 



(57) 
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The first two are obtained by linearizing the 3-divergence of their corresponding vector equations 



(41) and (43). The third results directly from definition (p2p) via the laws governing commutations 
between time derivatives and spatial gradients of scalars and spacelike vectors (see eqns ( |103| ) and 
(|104| ) or eqn (|105|) in Appendix B). Finally, the last is a simple rearrang ement of ©.0 



By combining equations (|54])-(|5"6|), or by linearizing the 3-divergence of ([5l]), we obtain the 
following second order differential equation for the evolution of the spatial matter aggregations 




v 3 ' c 2 



= 2(3) V 2 A 




^V 2 B. (59) 

The rest of the variables evolve in accordance with the propagation formulae, 

A+ «^ B . (60) 
3c z (l + w) 3/xc z (l + w) 

and 

4 46 / c 2 \ 

3(1 + w) 3(1 + w) \c z I 

where the latter is obtained by substituting Z in ([56]) from (|54|).P| This is the system that governs 
the evolution of spatial matter aggregations in a perturbed FFRW universe that contains a single 
barotropic perfect fluid of infinite conductivity and is permeated by a weak cosmological magnetic 
field. 

The equations obtained here are significantly simpler and more transparent than their vector 
counterparts of section 4.3, especially as far as the role of the magnetic field is concerned. The field 

12 The 3-divergence of (^) provides the evolution formula of C = fiK^V'Cj, the scalar associated with spatial 
inhomogeneities in the curvature of the spacelike regions. Its form, 

C = 4 9 C ; 9g \ (3) V 2 A+ 2Q f H \ (3) V 2 B, (58) 

verifies that C is time-invariant on large scales irrespective of the model's rotational beha viour. Notice that one 
immediately recovers (|57]) from ( |58[ ) on using definition ( p3[ ) and the commutation law ( |lOE| ) in Appendix B. 

"Equations @ and @ suggest that when w = 0, as it is the case in the dust era for example, then B = 
4A/3(1 +w). So, during these periods spherically symmetric spatial variations in the energy density of the magnetic 
field grow (or decay) proportionally to those in the energy density of the matter. 
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no longer exerts its influence through some complicated combinations of curl's and vector products, 
but simply via the spatial gradients of its energy density. Moreover, these are exactly the quantities 
that matter for structure formation purposes. 



6 Particular Solutions 

In [pi] we considered the evolution of density inhomogeneities during the post-equilibrium era, when 
the universe is filled with a non-relativistic perfect fluid (i.e. p = => w, c 2 s = 0). There, based 
on the nature of the evolution equation for A, we argued that the magnetic effects on the growth 
of large-scale material aggregations are relatively unimportant. Here, the existence of formulae 



and (61) will enable us to confirm, refine and extend these conclusions as well as to study the 



behaviour of the density contrast in the radiation era. 
6.1 Harmonic Analysis 



Following [16 1, |17]] and ||18|| , we harmonically decompose the inhomogeneity variable A by writing 



it in the form of the sum 

A = J2^ {n) Q {n \ (62) 

n 

with ^VjAW = 0, = and ^V 2 Q^ = -n 2 Q^/S 2 . Similarly, K, and B may be written as 

JC = J2lC {n) Q in \ and B = J2B (n) Q in \ (63) 

n n 

where ^Vj/C = ^V{B = 0. Notice that the harmonic eigenvalue (n) coincides with the comoving 
wavenumber (v) because of the spatial flatness of the background universe.^ 

Substituting results (62) and ( |63| ) into equations (|59|)-(|6l|) the harmonics decouple to provide 
the following autonomous system 

AM = -^ + ^-2^eAW + 




5 + f) A )^ W 

,2 /i o„2~ 



+ i + f A W> (64) 



3c z (l + w) 3c z (l + w) 



14 If the unperturbed universe has open spatial sections (i.e. k = —1) then n 2 — u 2 + 1, with v 2 > 0. Conversely, 
when the background model is spatially closed (i.e. k = +1) the associated relation is n 2 = v(y + 2), where now 



and the fundamental mode corresponds to v = 1, pq ], Jl9| ]. 
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and 

= ^rAM + f 4 " AH, (66) 

3(1 + w) 3(1 + w) yc z J 

where c\ = H 2 / \x is the Alfven speed characterizing the propagation of hydromagnetic waves. 
6.2 The Radiation Era 

When radiation dominates w = c 2 /c 2 = 1/3 and the energy density of the matter falls as [i = 
(see eqn fll3])), suggesting, together with equation (|24]), that the the Alfven velocity remains 
constant along the fluid-flow lines (i.e. c\ = 0). Ignoring the cosmological constant (i.e. A = 0), it 
is preferable to express equations (|6^)-(|66|) with respect to the scale factor, S(t), 



-2(l- J™-) AH + + *™) BM (67) 

5 dS 2 I 2kM r c 2 ) + c 2 c 2 [ L + 2KM R c 2 r ' 



d5 d5 



(69) 



In the long-wavelength limit (i.e. z/ — ► 0, or equivalently v 2 S 2 / kM r c 2 <C 1)P1 the above system 
reduces to 

S 2 ^- = 2AM + %/cH - %?W (70) 

5^ = aM + ^bM, (71) 

dBM _ dAM 

" dS ■ {U > 

and accepts a power-law solution of the form 

AW(5) = ^AHzS z , (73) 



15 During the radiation epoch the scale factor evolves as S = [it 1 / 2 , with /3 = (4kM_rc 4 /3) 1//4 . Considering a physical 
scale much larger than the horizon (i.e. \ p h y s S> dff), and taking into account that X p hys ~ .SAcom, A com ~ and 
d_H ~ ci, we find that v 2 S 2 / kWIrc 2 < 1 on large scales. Clearly, subhorizon scales are characterized by the reverse 
inequality. 



18 



where A^ are arbitrary positive constants. The parameter z satisfies the cubic equation 



z 2 + 



2c\ 



1 + 



3c| 
2c 2 



(74) 



which has three real roots provided that c\jc < 3/11, pl}| , given in trigonometric form by, [22], 



1 + 2V7 f 1 



3^/ 
14c 2 



cos 



9 + 2kir 



(75) 



with k = 0, 1, 2, and 



cos# 



10 



1 + ^4 



7V7 



1 



94 

4c 2 " 



(76) 



In the absence of a magnetic field (i.e. c\ = 0), expressions (|75|), (f76|) provide the standard solutions 
z = 0,2, —1 associated with a magnetic- free universe (see for example |23| or fT7|). However, the 
coupling between the field and the 3-curvature obscures the overall magnetic effect upon the growth 
of the density contrast. As ( |76| ) reveals the field increases the cosine term in ( |75"D but at the same 
time decreases this term's coefficient, with the net effect depending on the field's relative strength. 
To clarify the magnetic impact we consider the case of a spatially flat (i.e. K^ v ' = 0) perturbed 
universe. Then the system (|70|)-(|72]) reduces to 



,d 2 A^ 
dS 2 



(77) 



"dS" " "dS"' (78) 

while the density contrast evolves as 

AW(5) = A^V 1 + A^V 2 , (79) 
where A^, A^ are constants and 



When the field is absent we recover the familiar evolution law (i.e. Z\ = 2, Z2 = —1) of a magnetic- 
free universe. Generally however, the large-scale magnetic effect is to reduce the growth rate of the 
density contrast in proportion to its relative strength. 

Conversely, the evolution of short-wavelength (i.e. v — ► 00, or equivalently v 2 S 2 / kMlrc 2 S> 1) 
density aggregations is governed by the following set of equations 

j2a(» ,,2o2 o„2 q„2„2 c<2 

g 2 - 17 g A(") I ^rM 6v c A b B (u) 

dS 2 ~ kM r c 2 c 2 2kM/jc 4 ' 1 j 



19 



(82) 



"dS" " ~dS~' (83) 

The lack of a general analytic solution forces us to ignore the effects of the spatial curvature. In 
this case the remaining equations accept the solution 

^(S) = A« am [ ^=\[^M\ + Af> cos | + ^ I . (84 ) 





So, small-scale matter aggregations oscillate with period 2ttc\/kM^/v^J 1 + 3c\/2c 2 . Relative to 
the non- magnetized case (see [^3|] for example) , the excess pressure supplied by the field has simply 
increased the oscillation frequency of the density contrast. 

Conclusively, the presence of a cosmological magnetic field during the radiation era does not 
cause significant changes in the evolutionary patterns of the density gradients. However, as far as 
their actual growth is concerned, the field impact is evidently negative, although still secondary to 
the effects induced by the pressure of the relativistic matter, which still dominates their evolution. 

6.3 The Dust Era 

After the radiation era ends, dust dominates and the energy density evolves as \i = with 
Md = 0. Thus, for vanishing cosmological constant the scale factor changes as S = at 2 / 3 , where 
t measures the observer's time and a = (3kM£>c 4 /4) 1 / 3 . Also, = 2/t and // = 4/3/vC 4 t 2 . During 
this period the Alfven velocity falls as c\ = E/ai 2 / 3 , with E = 0, reflecting the fact that the 
magnetic energy density drops faster than that of the matter. So, relative to a reference frame 



comoving with the expanding fluid, equations (64)-(|6q) become 



dt 2 3t dt + 3t 2 + 9c 2 ai 8 / 3 3c 2 ai 8 / 3 { 4a 2 J ' 1 ' 

d£M 4E 



dt 3c 2 at 5 / 3 

dBM _ 4 dAM 
dt ~ 3 dt 



B^\ (86) 



(87) 



On superhorizon scales (i.e. v — > 0, or equivalently z^ 2 c 2 t 2 / 3 /a 2 -C 1)0 the above system 
simplifies into 

^ = _2_ 8E (y) _ 2E M 

dt 2 3t dt 3t 2 9c 2 at 8 / 3 3c 2 crf 8 / 3 

16 The post-equilibrium evolution of the scale factor implies that the long-wavelength condition \ p h ys 3> du trans- 
lates into u 2 c 2 t 2/s /a 2 < 1. 
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dt 



4E 



3c 2 at 5 /3 



fit 



3 (it 



(90) 



To obtain an analytic solution, we assume that the perturbed universe has flat spatial sections 
(i.e. K. = 0) and also consider the special case where fiM = 4AM /3. Then, we are left with the 
following differential equationp] 



d 2 A^ 



4 dAM 
3i di 



3t 2 



1 



4E 



3c 2 Crf 2 / 3 



A^. 



(91) 



Notice that at later times (i.e. i — * oo) the magnetic term in the parenthesis becomes completely 
irrelevant. So, in agreement with ffl], we recover the power-law evolution 



A W = A H r i + A M f 2/ 3j 



(92) 



also familiar from the study of a non-magnetized cosmological model. The alternative early time 
solution 



AM(t) 



» 



sm 



tV3 



+ A^cos / ' 



tV3 



et 1 / 3 + 



t 2/3 



(93) 



where e = 2\/2E/cy / a, suggests that under the magnetic influence the long-wavelength aggregations 
of the material component oscillate with an amplitude that increases as 

On scales well below the horizon (i.e. v — > oo, or equivalently i/ 2 c 2 t 2 / 3 /a 2 S> 1), equations 
-fl87|) become 



d 2 A (^) ^ 4 dA („) 2 8E 
~dW ~ ~3t~dT + 3^ + 9c 2 crf 8 / 3 



i/ 2 E 
2o^ 2 



B^\ 



(94) 



dl 



4E 



3c 2 at 5 /3 



(95) 



dBM _ 4 dAM 
cit ~~ 3 dt 



(96) 



17 According to equation (|9C|), the condition g^) /A^ = 4/3 requires that the same ratio holds at the initial 
moment too. In other words, solutions ( p2] ) and ( ]93| ) presume that, as the large-scale spatial variations in the 
magnetic and the fluid energy densities enter the post-equilibrium era, their ratio equals 4/3. Such a simplifying step 
is not unreasonable at all since, as equation (p9|) suggests, B^"' ~ A'"' by the end of the radiation era. 
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Again, by ignoring any effects from the spatial curvature we obtain the following power-law evolu- 
tion for the density contrast 

AM (t) = A { "k Zl + A^V 2 , (97) 

with, 

1 / / 24^ 2 e\ , , 

^— e^T-wJ- (98) 

Notice that in the absence of the magnetic field (i.e. when E = 0) we are left with the well known 
solution (i.e. z%2 = 2/3, —1) of the magnetic-free case. In quantitative agreement with Ruzmaikina 
and Ruzmaikin we find that the field presence reduces the growth rate of the inhomogeneities 
proportionally to the ratio t 2 ^ 3 H 2 / [i<? . 

We conclude by arguing that the presence of a cosmological magnetic field always opposes the 
growth of matter aggregations, either by forcing them to oscillate or by reducing their growth rate. 
The magnetic influence ceases only at the later stages of the dust era, when the relative strength of 
the field becomes negligibly small. It should be emphasized that result (|97| ) refers to wavelengths 
that lie within the horizon but are much larger than the Jeans length at the time. Otherwise the 
pressure effects of the ordinary non-relativistic matter become important preventing the density 
gradients from growing. This fact, together with the oscillatory nature of solution (|93]), suggests 
that earlier in the dust era any actual growth is confined to scales comparable to the horizon size 
at the time. 



7 Conclusions 

We have explored the influence of a primordial magnetic field upon the kinematical and the dy- 
namical evolution of perturbed cosmological models containing perfect fluids with non-vanishing 
pressure. We employed the Ellis-Bruni covariant and gauge-invariant formalism, first applied to the 
analysis of magnetized cosmologies in Q , to derive the full set of equations determining the linear 
evolution of an almost-FFRW universe containing a perfectly conducting medium. Relative to the 
dust era examined in BJ , the principal new complexities are due to changes in the observer's motion 
under the simultaneous action of the perturbed medium and the magnetic field. These changes are 



best seen in the different form of the momentum density conservation law (see eqn (14)), which in 
turn implies a modified acceleration for the fluid. In fact, this is the reason for essentially all the 
extra complications in the evolutionary patterns of the pre-equilibrium era. We have quantified the 
magnetohydrodynamical effects upon the kinematics and the dynamics of a universe dominated 
by a barotropic perfect fluid. We found an acceleration that depends on density gradients as well 
as on the gradients of the field. It is no longer normal to the field vector and can have subtle 
effects upon the evolution of fundamental cosmological parameters. Of particular interest is the 
first order magneto- geometrical contribution to the deceleration parameter. We show that, unlike 
ordinary matter which always slows the expansion down, the magnetic field can act as a driving 
force through its interaction with the geometry of the spatial sections. An analogous effect is found 
upon the Ricci scalar of the observer's instantaneous rest space. On large scales, the magnetic field 
tends to smooth out the curvature of the spacelike surfaces and restore their initial flatness. 
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As in Q, we were primarily interested in studying the growth of density inhomogeneities in a 
magnetized environment. Here, we have defined four scalar variables that measure spatial variations 
in the energy densities of the medium (A) and the magnetic filed (B), spatial inhomogeneities in 
the expansion (2), and deviations from the spatial flatness of the background universe (/C). We 
provide a system of four linear first order differential equations that describes the evolution of 
these disturbances and ultimately dictates the behaviour of spatial matter aggregations. We have 
obtained analytic solutions both at the long and at the short-wavelength limit during the radiation 
and the dust eras. In [jl| we argued for the relative unimportance of the field during the dust era 
and on scales that exceed the horizon at the time. Here, we were able to confirm and also refine 
those results. More specifically, we have found that any magnetic effects upon long-wavelength 
matter aggregations cease completely as the dust era enters its later stages. During this period 
the inhomogeneities grow exactly as those in a non-magnetized universe. Soon after equilibrium 
however, the extra pressure of the field could have forced the density gradients to oscillate, thus 
preventing them from growing. Nevertheless, the weakness of the field means that such large-scale 
oscillations are short-lived and that the epoch of unimpeded growth begins almost immediately 
after equilibrium. On scales smaller than the horizon, but larger than the associated at the time 
Jeans length, the disturbances undergo a power-law growth but at a slower pace relative to the 
magnetic-free case. The field pressure also affects the evolution of the density contrast during the 
radiation era. Here, it adds to the pressure of the relativistic matter and impedes any further 
gravitational clumping of the medium. On large scales we have found that the field inhibits the 
growth of the inhomogeneities by an amount proportional to its relative strength, whereas on 
subhorizon regions it increases the frequency of their oscillations. In the radiation era the magnetic 
effect supplements that from the pressure of the relativistic matter. During this period, the fate of 
small-scale inhomogeneities is affected by plasma processes || . 
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Appendices 

A Gauge-invariance for the 3-gradients of Spatial Vectors 

In (see Appendix A there) it was stated that the metric of an exact FRW or Bianchi-I spacetime 
can always be brought in the diagonal form 

gij = diag^oo, #11,522,533), (99) 

with its components being functions of proper-time only (i.e. gu = gu(t) - no summation over 
i). Based on this we then proceeded to prove the gauge-invariance of A4ij, the spatial tensor 
that describes the variations of the magnetic field vector as seen by two neighbouring fundamental 
observers. However, though our initial statement is correct within a Bianchi-I and a spatially flat 
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FRW cosmology (see for example ]^4j for verification), it can not be extended to spatially curved 
FRW spacetimes. So, the gauge-independence of the magnetic 3-gradients has been established 
simply within perturbed FFRW and Bianchi-I universes. What we actually showed in p] was 
that the spatial flatness of the aforementioned spacetimes (together with their zero rotation) is 
sufficient for the 3-gradients of any homogeneous spacelike vector field to be gauge-invariant. Next 
we argue that within the limits of a FRW model such a requirement is also necessary. Indeed, in 
the observer's rest space the Ricci identity takes the form (see Appendix B) 

^V^V^ = —Uijhfv, + l {3) R qkji v g , (100) 

where ViU 1 = 0. The above equation, which is presented here in its exact form, clearly states 
that in a non-rotating spacetime (i.e. LVij = 0) the 3-gradients of any spacelike vector vanish (i.e. 
^WiVj = 0) only when ^RijkqV 1 = 0. Contracting the latter over the indices j and q we obtain 
the new restriction ^R^v^ = 0. In a FRW spacetime ^Rij = ^'Rhij/S, which means that the 
gauge-invariance of ^VjUj requires ^R to vanish. This in turn ensures that ^Rijkq = and 
therefore the spatial flatness of the model. Clearly, the introduction of 3-vectors into the spatially 
isotropic Friedmannian cosmologies is only an approximation. Nevertheless, the 3-gradients of such 
a homogeneous spacelike vector, cannot be treated as gauge-independent variables unless the FRW 
cosmology is spatially flat. 



B Auxiliary Relations 

Following [0] we point out that generally the operator ^ 3 ^Vj cannot be treated as the standard 
covariant derivative of a 3-dimensional hypersurface because in a rotating spacetime the defect 
tensor does not vanish. Thus one cannot assume the usual commutation relations but should use 
expressions that include possible rotational terms. A selection of such formulae can be found in 



[14] (see Appendix A there). Here we present only those essential to our analysis. 

Commutations between the spatial gradients of scalars and spacelike vectors are respectively 
given by 

(^V^V j] f = -^u ij f, (101) 

and 

^V^V^fc = ~U3ijh§V q + - {3) RqkjiV q , (102) 

where / can be any scalar and V{ is a spatial vector (i.e. ViU 1 = 0). Equation (|102| ) is also regarded 
as the general expression of the 3-Ricci identity. Commutations between the spatial gradients and 
the time derivatives of these quantities are governed by 

{3) V,/ - V ( (3) V,/)' = ~f ai + ±8<»> V,/ + WV,/ (a\ + J) , (103) 

and 

(3, V^, - hfh? ( {3) VkV q )' = ^Q^ViVj, (104) 
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where the latter appears here in its linearized form and applies only to first-order (i.e. v $ = in 



the background) spacelike vectors. Commutator ( 104 ) provides an additional first order relation, 
which plays an important role in our analysis. In particular, assuming that ^Vj/ vanishes in the 
background, we may linearize the 3-divergence Of ( |103[) ) to obtain 



((3) V 2 /)' - (3) V 2 / = IjU - |( 3 'V 2 /, (105) 

recalling that A = ^V^a* to first order. The above is used to derive the evolution formula of B, 
the scalar that describes spherically symmetric changes in the energy density of the magnetic field. 
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